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Abstract 

We study migration of DNA molecules through a microchannel with a series of electric traps 
controlled by an ac electric field. We describe the motion of DNA based on Brownian dynamics 
simulations of a beads-spring chain. The DNA chain captured by electric field escapes due to 
thermal fluctuation. We find that the mobility of the DNA chain would depend on the chain 
length; the mobility sharply increases when the length of the chain exceeds a critical value, which 
is strongly affected by the amplitude of the applied ac field. Thus we can adjust the length 
regime, in which the microchannel well separates DNA molecules, without changing structure of 
the channel. We also present a theoretical insight into the relation between the critical chain length 
and the strength of binding electric field. 
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I. INTRODUCTION 



Gel electrophoresis is commonly used in biological science to separate out DNA molecules 
according to their size. A longer DNA has lower electrophoretic mobility in the random 
environment of a gel, the separation therefore would be possible using dc electric field. 
However, it is well known that gel electrophoresis loses its efficiency for DNA molecules 
longer than about 20k base pairs (bp), beyond which the length dependence of the mobility 
gradually disappears. Pulsed-field electrophoresis has been widely used for separation of 
longer DNA molecules. This method also has a limit, both with respect to speed and the 
size, with an upper limit about 10 Mbp. 

On the other hand, separation devices are required to incorporate into a mili-meter size 
chip with the developments of highly integrated bioanalvsis systems — the so called micro 
total analysis systems (/x-TAS) or lab-on-a-chip devices However introducing gels into 

microchannels may be a « problem Gel^ee separate technic gflflHQ there- 
fore have been proposed recently and attracted considerable amount of attention. Ajdari and 
Prost proposed a gel-free separation technique which uses free-flow electrophoresis together 
with a series of trapping by induced-dipole forces Q|. They argued from their theoretical 
analysis that, in their model, the mobility of a DNA chain monotonically decreases with 
its size. More recently, Han and Craighead have been investigated a designed microchannel 
(entropic trap arrays) which has alternating thin and thick regions (4 , 0] . They found that 
longer DNA molecules trapped in the thick regions for shorter time. Thus a shorter DNA 
has lower mobility. This fact has been confirmed theoretically j^] and numerically jjj. We 
note that, as mentioned in Ref. jj], the free energy landscape of the above two systems are 
quite similar, the relation between mobility and DNA size shows the clear contrast. This 
discrepancy may come from the fact that, to consider escape of DNA chains from a elec- 
trode, Ajdari and Prost regarded the chain as a point particle. Contrary, it is experimentally 
observed that when a DNA molecule escapes from an entropic trap, the deformation of DNA 
plays an important role J|. 

In this paper, we revisit the Ajdari and Prost's model and investigate the mobility of 
DNA chains along the channel, taking the deformation of DNA molecule into account. Our 
simulation verifies that longer DNA is trapped for shorter time and the mobility suddenly 
increases when the chain exceeds a critical length. Furthermore, we found that this critical 



2 



chain length is strongly affected by the voltage applied across the electrodes. Thus, we 
could adjust the regime in which the device well separates DNA chains without changing 
the structure of the channel. These numerical insights will help us to design practical 
separation devices. 

This paper is organized as follows. In Section II, our separation device and the simulation 
model are descriced. In Sec. III. A, result of our simulation is discussed, and in Sec. III.B, a 
scaling analysis is presented and compared with the simulation. In section IV, we summarize 
our results. 



II. THE SIMULATION MODEL 



To our simulation, we adopt a bead-spring model, which represents a DNA molecule as 
N beads or monomers of diameter a. Neighboring beads are connected by harmonic springs 
with constant k as 

where r is distance between neighboring beads, k-Q the Boltzmann constant and T the 
temperature. To take the excluded volume effect into account, all beads interact each other 
with a Weeks- Chandler- Andersen (WCA) potential 



V WCA {i 
knT 



(2) 

otherwise. 



where r is the inter-beads distance. 

Let us write the thermal fluctuation force acts on the ith bead as f]i(t). We suppose that 
the viscous drag is proportional to the bead velocity with friction constant (. A general 
theorem from statistical mechanics relates the random force and the friction constant as 

(f]i(t)r]j(t')) = 2k B T(5 ij I5(t - t'), (3) 

where / is the 3x3 identity tensor, 5ij the Kronecker delta and 8(t — t') the delta function. 
The random force in our simulation is set so as to satisfy this relation in a discretized manner. 

Schematic views in rcz-plane and x?/-plane of the channel are shown in Fig. ^ We simulate 
the motion of beads-spring chains inside this channel whose cross-section is a square with 
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side d. For x direction, periodic boundary condition with distance L = 9d is assumed both 
for the electric field and the chain motion unless mentioned particularly. We assume that 
the beads do not adhere to the channel's wall; the WCA potential also applied between the 
beads and the wall. In the corners, the two repulsive forces come from different directions 
are summed up. 

Each bead carries a electric charge g ne t- For strongly charged polyelectrolytes such as DNA 
molecules, the charge q net should be regarded as the renormalized one due to the counterion 
condensation. According to the Oosawa-Maining theory [lj]], the charge of DNA, immersed 
in monovalent salt, is neutralized by 76%. Hence we set, for the bead represents m base 
pairs, q nct ~ 2.0(1.0 — 0.76)me~ = 0.48me~ (e~ is elementary electric charge). 

Next we describe electric trapping of DNA molecules. A charged particle in solution is 
readily polarized by external electric field due to the cloud of counterions and the charge 
itself. In the presence of oscillating electric field E*, the polarization energy of this particle 
could be written as AU = a\E*\ 2 /2, with a as polarizability. This energy gives force 
/* = — grad(Af/), which drives the particles toward the regions of stronger field amplitude. 



Some experiments confirmed that the trapping o: 



DNA in the strongest field amplitude could 
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occur over the range of frequency 1Hz ~ 1kHz 

Ajdari and Prost theoretically considered the one- dimensional motion of a DNA molecule 
dragged by an static electric field. An ac electric field is applied perpendicular to the 
direction of DNA migration. In their model, it is assumed that the frequency is high enough 
and so that the molecule feel a ratchet like potential V(x) as depicted in Fig. 121 In our 
simulation, we set "parallel-plate condensers" along the channel as shown in Fig. Hand apply 
dc voltage across them. We also apply an uniform electric field E in x-direction, which 
drives the DNA chain along the channel. The potential <ft inside the channel is obtained 
by solving the Laplace equation V 2 = with the boundary element method under the 
boundary conditions of = ±AU on the electrodes and s • V0 = on the other wall, where 
s is the surface normal. The potential <ft obtained is approximately equivalent to the V(x) 
in x direction. The electric force acts on a bead is given as ff = q net {E — V0(Vj)} where 
7*j is position of i-th bead. Neglecting inertia, we write the equation of motion for i-th bead 
as 

Cn = fl* + ff + (4) 
where fl nt represents sum of the bead-bead and bead-wall interactions. 



The natural units of our simulation are defined as follows; we use a as the unit length, k-&T 
as the unit energy, £ = k^T / 'aq net as the unit electric field strength and r = (a 2 /ksT as the 
unit time. The channel size d = 20a. The spring constant k is set to 10 2 , in which condition 
the equilibrium length of the springs would be 0.85<r and no chain crossing occurs Q- We 
integrated the Eq. (J3J) with second order stochastic Runge-Kutta algorithms ^| using a 
time step At = 1.75 x 10~ 4 r. The random force is represented as rfc = {k^T/ a) a/ 6t/ Atipi, 
where ifri is a random vector which has independent components uniformly distributed on 
[—1,1]. We used Mersenne twister algorithms Q| to generate i/v Initial conformation of 
a chain is set by putting beads with interval a on a trajectory obtained from self-avoiding 
random walk. 

Finally, we mention how the parameters of our model can be related to those of an 



experiments. We first regard a = 50nm, order of the persistence length of DNA|15l.ll6|. then 
one bead corresponds to 167 base pairs. The electric field £ = 65V/cm under condition of 
T = 300K. Using the stokes' formula ( = 6irari with 77 as the viscosity of water, we obtain 
r = 5.1 x 10~ 4 sec. 

The present model neglects the effect of conterions and h ydr o dynamic interactions. Ad- 
equacy of the simplifications were precisely discussed in Ref.|l2|. This model has been suc- 
cessfully described the migration of DNA in gel ^0, 11] and in a channel with entropic traps 



We hence consider that the model may reproduce experimental results qualitatively 



and could fit to experimental results with a few adjustable parameters. 



III. DISCUSSION 



A. Results of simulation 



In this subsection, we show the results of our simulation for DNA migration in the 
microchannel. Here we set E = £x, where x is the unit vector parallel to the x axis, 
and describe the applied voltage across each electrode with a dimensionless parameter c = 
2AU/£d. The number of beads ranges from = 5 to = 10 3 , which corresponds to DNA 
molecules from 0.8 to 170 kbp. We focus on a dimensionless mobility jl = fi£a/r of a single 
isolated chain as a function of c and N. Figure 01 shows the snapshots of our simulation, in 
which a chain of A^ = 80 escapes from an electric trap under the conditions of c = 1.0 and 
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1.2. We observe that the chain is stretched during the migration to the next trap. Figure. 
H] shows trajectories for chains of length N = 20, 30, • ■ ■ , 160 at c = 1.0. The flat regions 
of the trajectories indicate the trapping of chains on a electrode. Shorter chains tend to be 
trapped for longer time and longer chains only slightly affected by the traps. To evaluate 
mobility, we fit the trajectory of which the chain has passed through a trap 10 times, by a 
linear function. This slope is the dimensionless mobility fi. 

We show the relation between fx and N in Fig. Each plot is the average of 10 trials. 
For c = 0.8, the mobility fx gently increases and for c > 1.0 a sharp increase of mobility 
appears. In the case of c > 1.2 short chains less than N ~ 10 2 are permanently trapped on a 
electrode and the resultant mobility would be zero. In the migration of long chains N ~ 800 
under the condition of c = 1.5 and 1.7, we observe that the mobility again rises. This is 
because the chain is stretched two or three times longer than the space between adjacent 
electrodes and the trapping by electrodes would be less efficient. In this regime, to obtain 
the theoretical understanding in the mobility is not easy. Thus this would be our future 
work. In the following subsection, we consider the dynamics of DNA chains trapped by only 
one electrode. 



B. Critical chain length iV c 

The microchannel could be efficient at separating DNA molecules in region where mobility 
fx sharply depends on chain length. Let us write the critical chain length as N c , at which the 
mobility shows the sudden increase. Figure El suggests that, by varying c, we can "regurate" 
the microchannel for the DNA moleucles we try to separate. In the following, we hence 
consider the relation between N c and c. 

Let us write a trapping lifetime of DNA chains as r t (c,N). The dimensionless mobility 



could be written as 



to + n 

where to is the transit time between two adjacent electrodes. Figure EJa) is a schematic 
view of a bead-spring chain trapped on a electrode, in which n < N beads of the chain 
have overcome a electric barrier AU. Each bead experiences the ratchet like potential V(x) 
as shown in Fig. El^b), where we set V(x)\ x= d = 0. To calculate trapping life time T t , we 
roughly estimate total free energy difference AF(n) due to the escape of n beads. Let r g 
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and l g be the typical sizes of a "blob" trapped on a electrodes. These two could be given by 



r g ~ 



a(N-n) Ur , (6) 



Z 9 ~ a{N-n) v K (7) 

Unfortunately, we cannot analytically calculate the exponents v r and v\. However, if we 
neglect the side walls at y — and d, our simulation evaluates z/ r ~ 0.44 and v\ ~ 0.37. 
Using the lengthes r 9 and l g , we can roughly estimate the electric potential energy of a blob 
trapped on a electrode £/biob as 

where A is constant. When n beads have overcome a barrier AU, the whole chain gains 

potential energy nq net AU and loses aq net £n 2 /2 as the n beads go down the slope. The 

potential energy of the chain therefore could be written as 

l g r 2 Q S n 2 
U(n) ~ Aq net — y — + nq net AU - oq nct £— , 
a 1 2 

= q nct a£ {A(N - nf + - y j , (9) 

where 7 = 2v T + 2^ ~ 1.25. Because the increase of entropic free energy is proportional to 
nksT, the total free energy difference AF(n) is given by 

AF(n) U{n) - U(0) 



knT kuT 



+ n 



cd 11 2 

= A{(N -n) 1 -N^} + —n + n. (10) 

2a 2 

In Fig. Et^b), we plot AF(n) as a function of ratio n/N. For c > 1 the AF(n) increases at a 
small n and reaches a peak at n = n*, while for c = 0.5 it monotonically decreases with n. 
Let AF max be the maximum of AF(n). Trapping occurs when d(AF)/dn\ n=Q > 0, and r t is 
given by 

r t ~ exp(AF max /A;Br). (11) 

We can calculate AF max by solving the following equation 

d AF(n) 



dn k B T 



n=n* 

-A lN ^(l- + g-„* + l = 0. (12) 



Here we assume that number of beads N is large enough, and neglect the order iV 7 2 and 
the higher terms. Then we obtain 



The sharp increase of jl appears when the free energy barrier AF max = 0. From FigJBfb), 
we know that this condition occurs only when n* = 0. We therefore obtain 



with A as a fitting parameter. To estimate the value of N c from simulation data, we adopted 
the following procedure: first, we obtained the maximum and minimum mobility (/i max , A^min) 
for a fixed c in region 5 < iV < 800. We write number of beads N c ' with this length the 
simulation gives mobility (yU m ax+A i min)/2. Figure |H1 shows that Eq. (fH^ is in good agreement 
with N c ' with one adjustable parameter A = 3.4. 

IV. CONCLUDING REMARKS 

We have presented Brownian dynamics simulations of DNA migration in channels with 
electric trappings. The chain length dependence of mobility observed in our simulation 
does not agree with the theory presented by Ajdari and Prost 3?]. We conclude that this 
discrepancy is due to their simplification, in which DNA molecules escapes from electric traps 
as a point particle, because our numerical simulations demonstrated that the deformation 
of DNA chain plays an important role when it escapes from the traps. 

Furthermore, the critical chain length N c is strongly affected by the voltage applied across 
the electrodes. From the rough consideration, we found a simple scaling relation between iV c 
and the parameters of the present system, which agrees with simulations with one adjustable 
parameters. 
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FIG. 1: A schematic views of a channel with electric traps in xy and yz-plane. The electrodes are 
square in shape and placed along the channel with interval 3d. 




FIG. 2: A sketch of the electric potential in Ajdari and Prost's channel. 



[a)c=1.0 (b>c=1.2 




FIG. 3: A DNA chain of N = 80 which escapes from a electric trap under the conditions of (a) 
c = 1.0 and (b) c = 1.2. The gray rectangular parts indicate electrodes. The insets in right hand 
side show a view in xy-plane. 
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FIG. 4: Trajectories of DNA chains of N = 20, 30, ■ ■ ■ , 160 (c = 1.0). We plot only the x-component 
of the center of mass. 
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FIG. 5: Mobility of DNA chains for c = 0.8, 1.0, • • • , 1.7. as the function of N. The errorbars 
show the maximum and minimum values of mobility in 10 traials. 
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FIG. 6: (a) A schematic view of a model DNA (bead-spring chain) escaping from a trap, (b) 
Electric potential energy experienced by a bead. 
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FIG. 7: The plot of Eq. fi® for c = 0.5, 1.0 and 1.5 (A = 3.4). 
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FIG. 8: The black circles represent the critical chain length estimated from our simulation. The 
solid line is a plot of Eq. Q14|) fitted to our simulation with parameter A = 3.4. 
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